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ABSTRACT 


Coastal  bathymetry  was  measured  using  wave  motion  as  observed  by  a  commercial 
satellite  imaging  system.  The  linear  finite  depth  dispersion  relation  for  surface  gravity 
waves  was  used  to  determine  nearshore  ocean  depth  from  successive  images  acquired  by 
the  WorldView-2  satellite  of  the  coastal  area  near  Camp  Pendleton,  California. 

Principal  component  transforms  were  performed  on  co-registered  images  and 
principal  component  four  was  found  to  very  effectively  highlight  wave  crests  in  the  surf 
zone.  Change  detection  images,  which  included  principal  component  four  from 
successive  images,  contained  both  spatial  and  temporal  information.  From  these  change 
detection  images,  wave  celerity  could  be  determined  and  depth  inversion  could  be 
performed. 

For  waves  farther  from  shore,  principal  component  four  no  longer  highlighted 
wave  crests.  Waves  could  be  resolved  within  a  single  RGB  composite  image  with 
equalization  enhancement.  The  wavelength  of  a  wave  above  a  known  depth  was 
measured  and  the  wave  period  method  was  used  to  determine  depth  for  other  waves  in 
the  propagation  direction  of  this  wave.  Our  depth  calculations  compared  favorably  to  our 
reference  bathymetry.  The  spatial  resolution  for  this  method  of  determining  depth  is 
higher  and  perhaps  more  accurate  than  our  reference  bathymetry,  particularly  in  the  surf 
zone. 
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I.  INTRODUCTION 

A.  PURPOSE  OF  RESEARCH 

Wave  and  tidal  processes,  and  severe  storms,  can  dramatically  change  nearshore 
bathymetry  over  a  short  time  period.  Nearshore  bathymetry  can  be  accurately  collected 
from  a  boat  using  sonar,  except  in  denied  waters  where  these  types  of  measurements 
would  not  be  allowed.  Figure  1  gives  an  illustration  of  the  sonar  bathymetry 
determination  method  (Gao,  2009): 


GPS  Antenna 


Figure  1.  Sonar  bathymetry  determination  method  (From  [Gao,  2009,  modified  from 

Tripathi  &  Rao,  2002]) 

Due  to  the  changeable  nature  of  the  nearshore  bathymetry,  the  usefulness  of  a 
sonar  survey  would  be  temporary.  Remote  sensing  measurements  from  the  air  may  not 
be  feasible  since  the  aircraft  could  be  in  danger  of  being  shot  down  prior  to  acquiring  data 
over  denied  territory.  Remote  sensing  from  Low  Earth  Orbit  (LEO)  can  be  used  over 
denied  territory  since  the  platform  is  safely  out  of  the  range  of  most  weapons  systems. 

A  space-based  electro-optical  payload  will  need  to  have  sufficient  resolution 
(~2m)  and  the  ability  to  take  multiple  images  of  the  same  coastal  location  in  short 
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succession  (6  to  12  s)  in  order  for  the  bathymetry  to  be  extracted.  This  researeh  foeuses 
on  estimating  water  depth  in  nearshore  environments  using  WorldView-2  multispeetral 
images  taken  in  rapid  suceession.  By  eomparing  the  same  wave  in  multiple  images  and 
knowing  the  exaet  time  of  image  acquisition,  water  depth  can  be  estimated  using  the 
linear  dispersion  relation  for  surfaee  gravity  waves. 

B.  SPECIFIC  OBJECTIVES 

The  objective  of  this  study  was  to  test  the  potential  of  determining  nearshore 
bathymetry  using  multiple  images  of  a  eoastal  area  taken  in  rapid  sueeession.  Being  able 
to  determine  bathymetry  in  this  manner  allows  for  aequiring  knowledge  of  denied  area 
nearshore  bathymetry  from  a  few  quick  images  from  a  LEO  spaceeraft.  Specifieally,  this 
work: 

1)  Uses  multiple  WorldView-2  multispeetral  images  (which  have  been 
registered  to  each  other)  taken  in  rapid  sueeession  of  the  coastal  area  near 
Camp  Pendleton,  California  to  identify  nearshore  waves 

2)  Uses  image  proeessing  techniques  in  the  Environment  for  Visualizing 
Images  (ENVl)  software  to  enhanee  wave  position  and  measure 
wavelengths 

3)  Calculates  depth  and  compiles  a  list  of  depth  data  points 

4)  Compares  depth  data  to  reference  bathymetric  contour  data  in  order  to 
determine  aeeuraey  of  estimates 

The  motivation  for  this  work  comes  from  a  desire  to  use  modem  remote  sensing 
techniques  to  determine  depths  without  using  water  transparency  methods.  Water 
transpareney  methods  suffer  from  absorption,  refleetion,  and  seattering  from  ehoppy  seas 
and  from  ocelusions  eaused  by  suspended  partieulates,  marine  plant  life,  and  dark  bottom 
sediments  (Williams,  1947).  The  approaeh  used  here  relies  only  on  being  able  to  resolve 
the  top  of  the  water  waves  and  uses  a  well-known  and  simple  relationship  between  wave 
celerity  and  depth.  Being  able  to  determine  water  depth  in  this  manner  has  applieation  to 
military  operations,  as  well  as  disaster  relief,  and  humanitarian  aid  efforts. 
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II.  BACKGROUND 


A.  LINEAR  FINITE  DEPTH  DISPERSION  RELATION  FOR  SURFACE 

GRAVITY  WAVES 

The  speed  of  shallow  water  waves  is  independent  of  wavelength  (^)  or  wave 
period  (T)  and  is  controlled  by  the  depth  of  water  (Knauss,  1996).  Waves  travel  faster  in 
deeper  water.  The  speed  of  deep  water  waves  is  independent  of  depth  and  is  determined 
by  X  and  T  (Knauss,  1996). 

Equation  1  shows  the  relationship  between  celerity  (surface  wave  phase  speed), 
wavelength,  and  water  depth: 

c"=^tanh(M)  (1) 

k 

Where  c  is  the  wave  celerity,  g  is  the  gravitational  acceleration  constant,  k  is  the 
wave  number  {k  =  27rl A),  and  h  is  the  water  depth.  This  equation  is  known  as  the  linear 
finite  depth  dispersion  relation  for  surface  gravity  waves.  There  are  approximations  that 
can  be  used  for  this  relation  that  depend  on  water  depth  and  wavelength  (Knauss,  1996). 
These  approximations  are; 

when  kh  <  0.33,  tanh(kh)  =  kh  (2) 

and  when  kh  >  1 .50,  tanh(kh)  =  1  (3) 

These  approximations  yield  dispersion  relation  equations  specific  for  shallow 
water  waves  and  deep  water  waves.  Equation  4  shows  the  equation  for  shallow  water 
waves  where  water  depth  is  dependent  on  wave  celerity.  Equation  5  shows  the  equation 
for  deep  water  waves  where  wave  celerity  is  dependent  on  wavelength. 

=  gh  (4) 


(5) 
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For  intermediate  wavelengths,  Equation  1  must  be  used.  These  equations  show 
that  the  wave  period  {T  =  T/c )  remains  constant  from  deep  water  to  shallow  water,  but 
wave  celerity  and  wavelength  decrease  proportionally.  Wave  height  (H)  also  changes 
from  deep  water  to  shallow  water,  and  will  further  complicate  the  dispersion  relation 
(Knauss,  1996). 

B.  HISTORY  OF  REMOTE  SENSING  DETERMINATION  OF  WATER 

DEPTH 

1.  Early  Investigations 

The  Normandy  invasion  by  allied  forces  in  World  War  II  presented  many 
logistical  challenges  including  how  to  land  hundreds  of  boats  on  beaches  controlled  by 
enemy  troops  (Williams,  1947).  Since  nearshore  coastal  areas  are  dynamic  and  can 
change  over  relatively  short  periods  of  time,  even  previously  surveyed  bathymetry  may 
not  be  good  enough  to  ensure  the  success  of  military  operations  where  lives  may  be  at 
risk.  In  W.W.  Williams’  seminal  paper  from  1947,  he  outlines  four  ways  of  attempting  to 
accurately  determine  bathymetry  remotely; 

a.  The  Waterline  Method 

Using  images  of  shorelines  taken  at  different  times,  contours  of  the 
waterlines  are  drawn.  The  highest  contour  will  correspond  to  high  tide,  while  the  lowest 
contour  will  correspond  to  low  tide.  In  ideal  conditions,  this  method  will  only  be  able  to 
produce  bathymetric  contours  for  the  difference  between  high  and  low  tides,  which  may 
only  correspond  to  a  few  feet.  Uncertainty  using  this  method  includes;  image  scaling  and 
rectification,  accurately  mapping  contours  to  tide  height  from  tables,  false  waterlines  due 
to  the  breaking  waves  flowing  up  on  the  beach,  and  problems  resolving  the  actual 
waterline  in  the  images  (Williams,  1947). 

The  waterline  method  was  used  extensively  during  World  War  II  and  was 
also  used  during  World  War  I  to  determine  beach  gradients  on  the  coast  of  Flanders 
(Williams,  1947).  The  utility  of  this  method  is  limited  since  only  very  shallow  depths 
can  be  mapped  where  tidal  processes  are  present.  For  any  body  of  water  that  experiences 

very  limited  tides  or  no  tides  at  all,  the  waterline  method  is  not  useful. 

4 


b.  The  Transparency  Method 

Incident  light  on  the  surface  of  the  water  will  be;  reflected,  absorbed, 
transmitted,  and  scattered  by  varying  amounts  (Olsen,  2007).  The  light  that  is  transmitted 
will  reflect  off  of  light  colored  sediments  and  be  transmitted  up  through  the  water  and 
toward  the  imaging  sensor,  in  very  shallow  water.  As  the  water  gets  deeper,  the  two-way 
attenuation  of  the  transmitted  light  energy  prevents  the  reflected  energy  from  the  bottom 
sediments  from  escaping  the  water.  When  this  happens,  it  is  no  longer  possible  to  see 
bottom. 

This  method  suffers  from  uncertainties  associated  with;  choppy  seas 
(increasing  scattering),  suspension  of  particulates  in  the  water  (preventing  light  from 
transmitting  through  entire  water  column),  and  dark  plant  and  sediment  materials 
(increasing  absorption  and  decreasing  reflection  of  incident  light).  This  method  was  used 
to  map  a  section  of  the  Seine  river  in  1943  (Williams,  1947). 

c.  The  Wave  Period  Method 

This  method  involves  studying  images  of  deeper  water  waves  and 
determining  their  wavelengths.  The  distance  from  shore  must  be  determined  in  order  to 
reference  known  bathymetry  to  determine  depth  where  the  wave  is  located.  Figure  2 
shows  a  plot  of  curves  used  to  determine  wave  period. 


A  -  Wavi  infect 

Figure  2.  Curves  relating  wave  period  to  depth  and  wavelength  (From  [Williams, 

1947]) 
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When  wavelength  and  depth  have  been  established,  a  plot  of 
mathematically  derived  curves  relating  wavelength  and  depth  to  period  can  be  referenced 
in  order  to  determine  wave  period.  Once  this  wave  period  has  been  determined,  different 
wavelength  waves  that  intersect  a  horizontal  line  drawn  on  the  plot  for  the  period  can  be 
directly  related  to  water  depth  (Williams,  1947). 

The  wave  period  method  requires  an  initial  known  depth.  Small 
inaccuracies  in  the  resultant  period  estimates  can  cause  large  errors  in  depth  calculation 
(Williams,  1947).  This  method  also  suffers  from  the  difficulty  in  accurately  measuring 
wavelengths  in  deeper  water  where  individual  waves  are  not  as  easily  distinguishable 
(Williams,  1947). 


d.  The  Wave  Velocity  Method 

The  celerity  of  an  ocean  wave  can  be  used  to  determine  depth.  In  order  to 
determine  celerity,  wavelength  must  be  measured  from  successive  images  with  a  known 
interval  between  them.  Figure  3  shows  a  plot  of  curves  used  to  determine  water  depth. 


Figure  3.  Curves  relating  water  depth  to  wavelength  and  celerity  (From  [Williams, 

1947]) 
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Deep  water  waves  travel  faster  than  shallow  water  waves  and  have  no 
dependence  on  water  depth.  As  waves  reach  depths  shorter  than  their  length,  they  slow 
down.  In  shallow  water,  water  depth  can  be  calculated  from  wave  celerity  and 
gravitational  acceleration  alone  (Williams,  1947). 

The  wave  velocity  method  is  useful  since  depth  can  be  calculated  from  a 
simple  mathematical  expression  once  celerity  is  known.  In  shallow  water,  waves  are  well 
defined  with  long  bright  crests  and  are  easier  to  distinguish  in  imagery.  With  multiple 
images  of  the  same  coastal  region  and  a  known  interval  between  images,  wave  celerity 
can  be  calculated  directly.  Cameras  with  clocks  were  not  always  available,  which 
decreased  the  usefulness  of  this  method  (Williams,  1947). 

The  wave  methods  described  above  suffer  from  uncertainties  including; 
accuracy  of  image  acquisition  times,  scale  and  registration  of  images,  complex  wave 
patterns  skewing  measurements,  and  low  image  resolution  (Williams,  1947).  Modern 
remote  sensing  systems  and  tools  can  be  used  to  overcome  some  of  these  limitations. 
Remote  sensing  spacecraft  typically  have  very  accurate  clocks,  and  image  acquisition 
times  are  embedded  in  image  metadata  files.  Images  can  be  registered  to  each  other 
using  software  tools  like  the  Environment  for  Visualizing  Images  (ENVI).  Commercial 
imagery  from  space  can  provide  up  to  50cm  resolution  improving  the  ability  to  resolve 
ocean  waves. 

2.  Modern  Studies 

Remotely  sensed  images  from  space  can  provide  a  nearly  comprehensive 
summary  of  wave  processes  for  relatively  large  areas  of  the  sea  surface  (Dalrymple, 
Kennedy,  Kirby,  &  Chen,  1998).  A  single  image  of  a  coastal  area  can  be  used  to 
determine  bathymetry  through  the  use  of  the  linear  dispersion  relation,  while  the  “lagged 
correlation  method”  where  multiple  images  taken  over  a  short  time  period  are  used 
together  can  be  more  useful  (Dalrymple  et  ah,  1998).  With  a  single  image,  wavelength 
must  be  measured  and  period  estimated.  Period  does  not  change  from  deep  to  shallow 
water,  but  it  must  be  estimated  for  the  deep  water  case.  With  multiple  images,  short  time 
samples  (~1  s)  are  needed  for  deep  water  depth  estimation. 
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Sequential  images  of  the  same  target  loeation  will  show  waves  of  various  wave 
numbers  and  wave  frequeneies  propagating  with  eharaeteristic  phase  speeds.  The  lag- 
correlation  method  uses  multiple  images  to  extract  the  wave  number  and  wave  frequency 
(Dalrymple  et  ah,  1998).  These  parameters  can  be  used  to  estimate  depth. 

The  linear  dispersion  relation  can  be  rearranged  to  express  depth  as  a  function  of 
X  and  y  position  (Dalrymple  et  ah,  1998): 


h{x,y) 


1 

k{x,y) 


tanh  ' 


(6) 


When  the  derivative  of  Equation  6  is  taken  (and  h(x,y)  is  divided  out),  the  error 
associated  with  depth  can  be  expressed  as  a  function  of  the  error  in  wave  number  (g(kh)) 
and  wave  frequency  (f(kh))  (Dalrymple  et  ah,  1998): 
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The  error  terms  f(kh)  and  g(kh)  double  any  error  in  wavelength  determination  in 
shallow  water,  and  increase  exponentially  with  kh  (Dalrymple  et  ah,  1998).  Figure  4 
shows  how  these  error  terms  increase  with  increasing  depth. 


to  diplh 


Figure  4.  Plot  of  error  functions  f(kh)  and  g(kh)  versus  kh  (From  [Dalrymple  et  ah, 

1998]) 
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There  is  more  error  assoeiated  with  caleulating  depth  in  deep  water  sinee  wave 
phase  speed  has  little  dependenee  on  depth.  For  deep  water,  any  error  assoeiated  with 
estimations  of  wave  phase  speed  or  wave  number  will  be  greatly  increased  (Dalrymple  et 
ah,  1998).  Another  source  of  uncertainty  comes  from  wave  heights  causing  depth 
overestimation. 

X-Band  marine  radar  images  can  also  be  used  to  estimate  water  depth  in 
nearshore  areas.  A  sequence  of  radar  images  can  be  analyzed  to  map  how  wave  behavior 
changes  in  these  areas.  Wave  period  must  be  calculated  from  the  radar  data  and  used  to 
determine  water  depth  (Bell,  1999).  Wave  celerity  and  direction  can  be  mapped  with  the 
radar.  Bathymetric  inversion  is  used  once  peak  wave  period  has  been  measured  from  the 
radar  data.  Linear  wave  theory  use  for  nearshore  areas  suffers  from  uncertainties 
associated  with  nonlinear  processes,  but  allows  for  a  good  approximation  of  depth  (Bell, 
1999).  Figure  5  shows  a  plot  of  wave  celerity  using  an  X-Band  marine  radar  system. 
The  water  depth  plot  shown  in  Figure  6  was  calculated  using  the  radar  data  and  the  linear 
wave  theory. 
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Figure  5.  Plot  of  wave  eelerity  from  X-Band  marine  radar  (From  [Bell,  1999]) 
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Figure  6.  Plot  of  water  depth  calculated  from  X-Band  radar  data  (From  [Bell,  1999]) 


In  nearshore  areas,  the  wave  period  has  less  influence  on  wave  celerity. 
Nonlinear  wave  processes  will  contribute  to  larger  wave  celerities  than  predicted  with  the 
linear  dispersion  relation  (Bell,  1999).  This  will  cause  depth  estimates  to  be  greater  than 
actual  depths.  In  order  to  improve  this  uncertainty,  the  use  of  higher  order  wave  theory 
must  be  used  (Bell,  1999). 

One  known  previous  attempt  has  been  made  to  measure  coastal  bathymetry  using 
satellite  imagery  of  wave  motions  (Abileah,  2006).  This  implementation  uses  one  meter 
imagery  from  the  IKONOS  satellite  of  the  Coronado  Island,  San  Diego  area.  The  Fourier 
transform  is  used  to  convert  image  intensity  into  wave  number  (Abileah,  2006).  The 
wave  number  spectrum  at  a  time  t  is  related  to  the  wave  number  spectrum  at  time  t=0  by 
the  following  equation  where  U  is  ocean  current: 

-‘(jsN tanh([t|c/)+[c/, ,Uy ]  ) 

S^=Sge  ^  ^  (8) 
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The  phase  ehange  between  two  or  more  of  the  image  transformations  can  be  used 
to  determine  both  current  and  depth.  Figure  7  shows  an  example  of  how  this  process  can 
produce  depths  using  the  IKONOS  data  (Abileah,  2006): 


Figure  7.  Bathymetry  using  IKONOS  imagery  (left)  (color  scale  is  red-blue=0-22m 
depth)  and  corresponding  USGS  topographical  map  (right)  (From  [Abileah, 

2006]) 


Video  image  processing  can  also  be  used  to  estimate  the  depth  of  nearshore  areas. 
Amplitude  and  phase  are  computed  from  the  video  imagery  and  the  slope  of  the  phase  is 
used  to  calculate  the  cross-shore  wave  number  (Stockdon  &  Holman,  2000).  Depth  can 
then  be  calculated  using  a  modified  depth  inversion  equation  where  kx(x)  is  cross-shore 
wave  number,  ky(x)  is  long-shore  wave  number,  and  is  variance  (Stockdon  &  Holman, 
2000): 


h{x)  = 


tanh  ' 

( 
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(7^ 

\esxf*e,(xf 

(9) 


Pixel  intensity  as  a  function  of  position  and  time  is  used  as  the  input  signal  for  the 
video  technique.  The  data  is  Fourier  transformed,  normalized  by  a  matrix  of  variances  at 
each  position,  and  the  complex  eigenvector  of  the  first  mode  is  used  for  amplitude  and 
phase  calculations  (Stockdon  &  Holman,  2000).  Figure  8  shows  sample  results  from  this 
process  (Stockdon  &  Holman,  2000): 


11 


Figure  8.  Plot  of  amplitude,  phase,  cross-shore  wave  number,  and  depth  calculated  from 

video  data  (From  [Stockdon  &  Holman,  2000]) 


Solitary  theory  modifies  the  linear  dispersion  theory  to  account  for  wave  heights. 
If  wave  heights  can  be  estimated,  the  celerity  of  a  solitary  wave  (single  wave  crest)  in 
shallow  water  is  just  (Stockdon  &  Holman,  2000): 

c  =  4g{h  +  H)  (10) 


The  depth  uncertainty  between  using  linear  theory  versus  solitary  theory  can  be 
shown  by  relating  the  two  equations  to  one  another  where  Cs  is  wave  celerity  from 
solitary  theory,  and  H  is  wave  height  where  H/h  is  assumed  to  be  0.42  in  the  surf  zone 
(Stockdon  &  Holman,  2000): 
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(11) 


This  equation  shows  that  errors  associated  with  depth  calculations  using  linear 

theory  of  up  to  42%  may  be  possible.  In  order  to  mitigate  these  errors  without  including 
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wave  height  (which  may  be  difficult  to  obtain  from  remotely  sensed  imagery)  in  depth 
calculations,  low  amplitude  waves  a  little  farther  away  from  shore  should  be  used 
(Stockdon  &  Holman,  2000). 

The  linear  dispersion  relation  has  been  found  to  very  accurately  predict  water 
depth  for  waves  outside  of  the  surf  zone  using  extensive  field  measurements  (Holland, 
2001).  Inside  the  surf  zone,  wave  height  increases  celerity.  Using  linear  theory,  which 
does  not  account  for  wave  height,  will  underestimate  the  celerity  and  result  in 
overestimates  for  water  depth  (Holland,  2001).  The  surf  zone  of  a  given  location  can 
vary  greatly.  The  surf  zone  at  a  site  near  Duck,  NC  was  found  to  be  where  water  was  less 
than  4m  deep  (Holland,  2001).  For  depth  inversion  calculations  outside  the  surf  zone, 
errors  of  3%  to  9%  from  measured  depths  were  observed.  Inside  the  surf  zone,  errors 
over  50%  were  common  (Holland,  2001).  Figure  9  shows  how  errors  increased  with 
shallower  depths: 
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Figure  9.  Plot  of  predicted  versus  measured  depth  for  multiple  experiments  (From 

[Holland,  2001]) 
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Taking  wave  height  into  account  when  performing  bathymetric  inversion  will 
improve  the  accuracy  of  shallow  water  depth  prediction.  The  problem  remains  to  be  able 
to  determine  wave  height  from  remote  sensing  imagery. 

If  wave  heights  cannot  be  determined  from  remote  sensing  imagery,  they  could  be 
estimated  in  order  to  reduce  uncertainty.  The  depth  of  water  would  be  that  calculated  by 
linear  dispersion  theory  minus  the  height  of  the  wave.  Concentrating  on  resolving  waves 
outside  of  the  surf  zone  allows  a  solution  that  neglects  wave  height,  but  will  have 
challenges  associated  with  being  able  to  resolve  the  waves  in  optical  imagery  (Holland, 
2001). 

Studies  have  been  performed  with  multiple  images  from  airborne  platforms  to 
determine  nearshore  ocean  depth  (Dugan,  Piotrowski,  &  Williams,  2001).  The  images 
were  mapped  and  rectified,  and  radiance  data  was  Fourier  transformed  to  create 
frequency-wave  number  spectra  for  the  modulations  caused  by  ocean  waves.  This 
produced  spectra  with  information  on  the  wave  characteristics  (celerity,  wavelengths, 
etc.)  that  were  used  to  predict  water  depth  using  linear  dispersion  theory  (Dugan  et  al., 
2001). 

In  order  to  overcome  problems  associated  with  image  registration,  an  inertial 
navigation  system  (INS)  was  mounted  on  the  camera  turret  in  order  to  provide  camera 
attitude  and  camera  position  using  the  Global  Positioning  System  (GPS)  (Dugan  et  al., 
2001).  The  collection  geometry  of  this  airborne  platform  is  shown  in  Figure  10  (Dugan 
et  al.,  2001): 
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Figure  10.  Airborne  nearshore  area  collection  geometry  a)  perspective  view  b)  overhead 

view  (From  [Dugan  et  al.,  2001]) 

The  process  of  producing  high  resolution  frequency-wave  number  spectra  was 
enabled  primarily  through  the  use  of  large  image  sub-patches,  high  quality  camera 
hardware,  and  accurate  mapping  techniques.  This  resulted  in  depth  estimates  with  errors 
of  5%  to  13%  with  respect  to  the  surveyed  bathymetry  (Dugan  et  al.,  2001). 

Uncertainty  in  the  depth  inversion  process  can  be  characterized  by;  measurement 
accuracy  for  wave  phase  speed  and  wave  height,  and  choice  of  depth  inversion  model  for 
the  observed  wave  conditions  (Catalan  &  Haller,  2008).  The  complex  relationship 
between  wave  phase  speed  and  water  depth  involves  nonlinear  processes,  which  can  be 
very  difficult  to  account  for  with  remotely  sensed  data  (Catalan  &  Haller,  2008).  Linear 
models  of  wave  celerity  have  been  found  to  underestimate  phase  speed  since  some  of  the 
wave  energy  is  going  into  the  z  direction  (wave  crest  and  trough)  (Catalan  &  Haller, 
2008).  The  two  main  things  that  affect  phase  speed  estimation  are;  amplitude  dispersion, 
which  accounts  for  nonzero  wave  amplitudes,  and  frequency  dispersion,  which  accounts 
for  the  relative  depth  of  water  (Catalan  &  Haller,  2008). 

Large  scale  laboratory  experiments  have  been  performed  using  high  resolution 
remote  sensing  video  data,  and  surface  elevation  data  from  in-situ  wave  gauges  in  order 
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to  compare  the  accuracy  of  wave  celerity  models  in  nearshore  areas  (Catalan  &  Haller, 
2008).  Figure  1 1  shows  example  video  data  from  these  experiments: 
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Figure  1 1 .  Example  video  data  showing  a)  intensity  timestack  and  b)  pixel  intensity  time 

series  (from  [Catalan  &  Haller,  2008]) 


The  surf  zone  was  chosen  for  comparison  since  in  this  area,  nonlinearities  are 
known  to  add  uncertainty  in  the  depth  inversion  method  of  extracting  bathymetry. 

- c 

Results  for  this  experiment  are  shown  in  Table  1  where  R  is  the  mean  relative  error  and 

— c 

Rrms  is  the  average  root-mean-square  error  for  each  model  (Catalan  &  Haller,  2008). 
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Table  1.  Experimentally  derived  errors  for  10  different  wave  celerity  models  (From 

[Catalan  &  Haller,  2008]) 


Model 

Equation 

^RMS 

Linear  theory 

=  1  tanh  {kh) 

-7.4 

10.5 

Solitary 

+ 

II 

17.9 

19.2 

Modified  shallow 

c=1.3v^ 

24.0 

25.7 

Cnoidal  theory 

/  r  .  /  -\  1 

7.0 

9.1 

c=^gh[l+^(2-m-3^)] 

Modified  cnoidal 

c  =  >/?*(!  +f{m)H'/h) 

5.0 

10.1 

13.3 

14.7 

9.1 

11.9 

Bore 

13.3 

14.7 

Shock 

c  =  -2v^  +  2y/^  +  ^(*1  (H^)) 

24.1 

25.1 

KD86 

=g/k(l  + /i6^D)tanh(A/i  +^€) 

2.8 

8.0 

Hedges 

=  1  tanh  {k{h  H-  ^)) 

11.1 

12.7 

Booij 

c2=|tanh(ib(/*  +  H^/2)) 

2.7 

8.0 

Several  of  the  nonlinear  models  from  the  academic  literature  (KD86  (Kirby  & 
Dalrymple,  1986)  and  Booij  in  particular)  perform  better  than  the  linear  models  in  the 
surf  zone  (Catalan  &  Haller,  2008).  The  standard  linear  dispersion  relation  for  surface 
gravity  waves  or  “Linear  theory”  is  shown  to  underestimate  phase  speed  by  not  taking 
wave  height  into  account.  The  “Modified  shallow”  model,  which  is  a  version  of  the 
“Solitary”  model  where  the  relation  H/h  is  taken  to  be  a  constant  of  0.42,  is  shown  to 
greatly  overestimate  the  wave  phase  speed.  All  of  the  models  used,  except  for  “Linear 
theory”  and  “Modified  shallow,”  require  a  known  wave  height  or  local  water  depths 
under  wave  crest  and  trough  in  order  to  obtain  more  accurate  results  (Catalan  &  Haller, 
2008).  The  determination  of  wave  heights  is  especially  difficult  to  derive  from  remote 
sensing  imagery  because  the  relative  height  of  the  waves  is  very  small  with  respect  to  the 
platform  to  target  range. 

The  linear  theory  can  be  compared  to  KD86  (the  most  complex  celerity  equation) 
in  terms  of  how  accurately  they  can  determine  depth.  Table  2  shows  experimental  data  of 
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depth  inversion  using  each  model  where  a  positive  percent  error  corresponds  to  an 
overprediction  of  depth,  D  refers  to  difference  errors  and  R  relative  errors  for  multiple 
cases  (Catalan  &  Haller,  2008): 


Table  2.  Experimental  error  estimates  of  linear  and  KD86  models  (From  [Catalan  & 

Haller,  2008]) 


Case 

Linear 

KD86 

D 

R 

D 

R 

m 

m 

% 

% 

m 

m 

% 

% 

35 

0.05 

0.10 

9.41 

19.26 

-0.10 

0.14 

-13.13 

20.47 

36 

0.05 

0.11 

5.28 

14.36 

-0.12 

0.16 

-20.29 

27.20 

37 

0.25 

0.26 

41.10 

44.83 

0.08 

0.11 

14.75 

19.85 

38 

0.23 

0.25 

39.02 

44.89 

0.05 

0.11 

10.93 

20.55 

39 

0.15 

0.18 

27.15 

33.30 

-0.03 

0.11 

-1.39 

15.16 

Median 

0.17 

0.18 

28.20 

31.65 

-0.01 

0.05 

0.94 

9.70 

Table  2  shows  experimentally  that  the  KD86  model  can  be  more  accurate  on 
average  at  predicting  depth  in  the  surf  zone.  The  greater  accuracy  comes  at  the  price  of 
complexity  and  the  need  for  knowing  the  wave  heights.  The  full  KD86  model  is 
described  by  the  following  equations  (Catalan  &  Haller,  2008): 

=  gj  A:(l  +  yjf^D^anh^M  +  (12) 

where  e  =  kHjl  and  H  is  the  wave  height  and: 


8  +  cosh -  2  tanh  kh 

8  sinh"'  kh 

(13) 

(kh^=  tanh^  (M) 

(14) 

rnu^  ^  kh  V 

(15) 

If  wave  heights  are  known,  a  wave  celerity  model  that  incorporates  this  factor  can 
be  used  in  the  depth  inversion  process  for  determining  water  depth.  In  all  other  cases,  a 
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linear  model  must  be  used.  Figure  12  shows  a  flowchart  for  the  depth  inversion 
algorithm  where  dashed  lines  are  the  steps  for  linear  inversions  and  solid  lines  are  the 
extra  steps  required  for  nonlinear  inversions  (Catalan  &  Haller,  2008): 
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Figure  12.  Depth  inversion  algorithm  flowchart  (From  [Catalan  &  Haller,  2008]) 
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III.  EXPERIMENTAL  APPROACH 


A.  PROBLEM  DEFINITION 

Depth  inversion  using  space  based  imagery  is  a  convenient  method  to  measure  the 
bathymetry  of  nearshore  areas  in  remote  or  denied  locations.  This  method  has 
advantages  over  the  transparency  method,  which  suffers  from  absorption,  reflection  and 
scattering  (Olsen,  2007)  off  of  suspended  particulate  matter,  aquatic  plants,  and  dark 
bottom  sediments  (Williams,  1947).  Depth  inversion  requires  that  there  are  visible  waves 
in  the  data,  and  that  their  wavelengths  and  celerities  can  be  determined. 

This  study  uses  multispectral  imagery  acquired  by  the  WorldView-2  spacecraft  of 
the  coastal  area  near  Camp  Pendleton,  California  to  see  if  it  offers  advantages  for 
determining  nearshore  depth  using  the  linear  finite  depth  dispersion  relation  for  surface 
gravity  waves.  Commercial  imaging  spacecraft  can  provide  up  to  50  cm  resolution 
panchromatic  imagery,  but  this  is  excessively  high  when  used  for  resolving  things  as 
large  as  ocean  waves.  Because  increased  spatial  resolution  will  not  improve  the  ability  to 
perform  linear  depth  inversion,  we  will  determine  if  the  increased  spectral  resolution 
provided  by  WorldView-2  can  provide  a  benefit  for  this  method. 

B.  MATERIALS 

1.  WorldView-2  Sensor 

WorldView-2  is  the  third  satellite  in  DigitalGlobe’s  commercial  satellite 
constellation  (Figure  13).  It  was  built  by  Ball  Aerospace  &  Technologies  Corp.  and 
launched  in  October  2009. 
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Figure  13.  WorldView-2  Satellite  (From  [DigitalGlobe,  2009]) 


WorldView-2  is  the  first  commercial  satellite  to  combine  both  high  resolution 
panchromatic  and  8-band  multispectral  sensor  capabilities  and  is  a  significant  spectral 
performance  improvement  over  DigitalGlobe’ s  two  previous  satellites  QuickBird  and 
Worldview- 1  (Figure  14).  The  panchromatic  sensor  is  capable  of  46  cm  resolution  at 
nadir.  The  multispectral  sensors  are  capable  of  1.84  m  resolution  at  nadir  (DigitalGlobe, 
2009). 
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Figure  14.  Comparison  of  DigitalGlobe  satellite  spectral  coverage  (From  (DigitalGlobe, 

2010b]) 


Table  3  shows  important  specifications  of  the  WorldView-2  spacecraft  and 
payload: 
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Table  3.  WorldView-2  Specifications  (From  [DigitalGlobe,  2009]) 


Launch  Information 

Date:  October  8, 2009 

Launch  Vehicle:  Delta  7920  (9  strap-ons) 

Launch  Site:  Vandenberg  Air  Force  Base,  California 

Orbit 

Altitude:  770  kilometers 

Type:  Sun  synchronous,  10:30  am  descending  node 

Period:  100  minutes 

Mission  Life 

7.25  years,  including  all  consumables  and  degradables  (e.g.  propellant) 

Spacecraft  Size, 

Mass  and  Power 

4.3  meters  (14  feet)  tall  x  2.5  meters  (8  feet)  across 

7.1  meters  (23  feet)  across  the  deployed  solar  arrays 

2800  kilograms  (6200  pounds) 

3.2  kW  solar  array,  100  Ahr  battery 

Sensor  Bands 

Panchromatic:  450  -  800  nm 

8  Multispectral: 

Coastal:  400- 450 nm  Red:  630-690  nm 

Blue:  450  -  510  nm  Red  Edge:  705  -745  nm 

Green:  510  -  sBonm  Near-IRi:  770-895  nm 

Yellow:  585  -  625nm  Near-IR2:  86o-io4onm 

Sensor  Resolution 

Panchromatic:  0.46  meters  GSD  at  nadir*,  0.52  meters  GSD  at  20°  off-nadir 
Multispectral:  1.84  meters  GSD  at  nadir*,  2.08  meters  GSD  at  20°  off-nadir 

Dynamic  Range 

11-bits  per  pixel 

Swath  Width 

16.4  kilometers  at  nadir 

WorldView-2  is  also  capable  of  collecting  multiple  point  targets  that  are  16.4  km 
X  16.4  km  in  size.  Multiple  images  of  the  same  point  target  can  be  collected  in  rapid 
succession  due  to  the  increased  agility  of  this  spacecraft  (acceleration=1.5  deg/s/s, 
rate=3.5  deg/s)  (DigitalGlobe,  2009). 

2.  The  Environment  for  Visualizing  Images  4.7  (ENVI) 

The  Environment  for  Visualizing  Images  version  4.7  was  used  to  manipulate 
WorldView-2  imagery  of  nearshore  waves.  ENVI  4.7  is  a  powerful  tool  used  for  image 
processing  for  remote  sensing  applications  (ITT  Visual  Information  Solutions,  2010).  It 
is  capable  of  performing  image  to  image  registration.  Image  features  can  be  accurately 
measured  using  the  built-in  measurement  tool.  Principal  components  transformations  can 
be  performed  on  images  for  decorrelation  of  variables.  A  form  of  change  detection  can 
be  performed  on  co-registered  images  by  representing  different  images  as  different  colors 
in  an  RGB  triple  representation  of  the  data. 
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C.  METHODS 

1.  Image  Registration 

The  World  View-2  imagery  provided  by  DigitalGlobe  of  the  coastal  area  near 
Camp  Pendleton  are  ortho-ready  standard  2A  format  and  are  in  geographic 
latitude/longitude  coordinates.  In  this  format,  the  image  is  projected  onto  a  reference 
ellipsoid  without  the  rough  terrain  model  applied.  A  constant  base  elevation  is  used 
instead  (DigitalGlobe,  2010a).  Image  to  image  registration  using  ground  reference  points 
can  be  performed  on  these  images  using  the  ENVI  map  registration  function.  Images  of 
the  Camp  Pendleton  area  were  co-registered  using  the  resampling,  scaling  and  translation 
(RST)  method  with  nearest  neighbor  resampling  for  ocean  wave  comparison  purposes. 

2.  Principal  Components  Transform 

Using  principal  components  analysis,  a  number  of  potentially  correlated  variables 
can  be  linearly  transformed  into  a  smaller  set  of  uncorrelated  variables  (Jackson,  2003). 
The  increased  degrees  of  freedom  offered  by  the  eight  band  multispectral  image  data  can 
be  examined  using  image  processing  techniques.  Principal  components  transforms  were 
performed  and  statistics  computed  on  the  multispectral  images  using  the  ENVI  principal 
components  transform  function  in  order  to  determine  whether  these  could  aid  in  wave 
detection. 

3.  Change  Detection 

The  WorldView-2  image  data  consist  of  multiple  time  spaced  images  of  the  same 
coastal  scene.  Co-registered  images  can  be  placed  in  each  band  of  an  RGB  triple  image 
representation  to  show  how  things  change  from  image  to  image.  For  comparing  two 
images,  the  second  image  is  placed  in  the  red  and  green  bands  while  the  first  image  is 
placed  in  the  blue  band.  Multiple  images  were  compared  in  this  way  to  show  the  wave 
moving  from  blue  to  yellow  (red  combined  with  green)  from  one  image  to  the  next 
image. 
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IV.  OBSERVATIONS 


A.  WORLDVIEW-2  IMAGERY  OF  CAMP  PENDLETON 

The  WorldView-2  imagery  of  the  coastal  area  near  Camp  Pendleton,  California 
provided  by  DigitalGlobe  consists  of  ten  panchromatic  and  multispectral  images  taken  in 
rapid  succession  during  a  single  pass  of  the  spacecraft  over  the  area  on  March  24,  2010. 
The  images  are  in  ortho-ready  standard  2A  format  and  are  in  geographic 
latitude/longitude  coordinates.  Figure  15  shows  how  the  images  are  broken  up  into  two 
parts  each  and  designated  at  row  1  column  1  (RlCl)  and  row  2  column  1  (R2C1). 


Figure  15.  DigitalGlobe  ortho-ready  standard  2 A  image  format  (From  [DigitalGlobe, 

2010a]) 

The  Camp  Pendleton  area  scene  consists  of  about  half  ocean  and  half  land.  This 
improves  the  ability  to  co-register  the  images  while  providing  a  wide  field  of  view  of  the 
ocean  to  find  different  types  of  waves.  Figure  16  shows  the  layout  of  the  Camp 
Pendleton  scene  with  both  row  1  and  row  2  in  a  Google  Earth  representation  of  the  first 
collected  image. 
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Figure  16.  Google  Earth  representation  of  WorldView-2  image  data  for  Camp  Pendleton 

area 

Image  metadata  for  the  ten  multispectral  images  were  examined  for  information 
on  image  acquisition  times.  The  “first  line  time”  corresponds  to  exposure  of  the  first  line 
in  the  image.  The  time  delta  between  images  is  shown  in  Table  4.  The  time  between 
successive  images  varies  from  10.40  seconds  to  11.00  seconds.  The  time  delta  between 
the  first  and  second  image,  second  and  third  image,  and  third  and  fourth  image  is  the 
same  to  a  hundredth  of  a  second.  The  times  given  in  Table  4  are  Zulu  time.  Pacific 
daylight  savings  time  began  on  March  14  in  2010,  so  1900  corresponds  to  noon  PDT. 
Note  that  the  image  (folder)  names  are  arbitrary,  and  derive  from  the  vendor’s  processing 
string. 
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Table  4.  Multispectral  image  directories  showing  time  between  images 


World  View-2  Data:  March  24,  2010 

Directory  Name:  052365928020 01 

Folder  Name 

Mean  GSD  X 

Mean  GSD  Y 

Mean  Sat  AZ 

Mean  Sat  EL 

First  Line  Time 

Delta  Time 

P002 MUL 

2.087 

2.374 

277.8 

61.2 

19:03:31.83 

0 

P007 MUL 

2.102 

2.425 

267.1 

60.1 

19:03:42.63 

10.80 

P009 MUL 

2.149 

2.503 

257.6 

58.2 

19:03:53.43 

10.80 

P006 MUL 

2.229 

2.609 

249.5 

55.8 

19:04:04.23 

10.80 

P008 MUL 

2.337 

2.741 

242.9 

53.1 

19:04:14.83 

10.60 

P004 MUL 

2.481 

2.910 

237.3 

50.1 

19:04:25.83 

11.00 

P010 MUL 

2.647 

3.107 

232.9 

47.2 

19:04:36.43 

10.60 

P005 MUL 

2.841 

3.338 

229.2 

44.4 

19:04:47.03 

10.60 

P003 MUL 

3.057 

3.603 

226.2 

41.6 

19:04:57.43 

10.40 

P001_MUL 

3.299 

3.908 

223.6 

39.1 

19:05:07.83 

10.40 

The  images  can  all  be  registered  to  either  image  eight  or  image  four  which  are  the 
middle  two  images.  The  first  three  images  were  registered  to  image  eight  using  four 
ground  reference  points  located  near  the  water  line.  The  first  ground  reference  point  was 
chosen  near  the  upper  left  comer  of  the  image,  the  second  was  chosen  near  the  bottom  of 
the  image,  and  the  third  and  fourth  were  chosen  near  the  center  of  the  coastline  in  the 
image.  Figure  17  shows  an  example  image  with  the  ground  control  points.  A  limited  set 
of  control  points  was  purposely  chosen  because  the  focus  was  on  the  flat  ocean  surface. 
It  is  difficult  to  estimate  the  registration  error  for  the  open  water,  given  the  lack  of  control 
points  in  the  water.  Typical  errors  for  this  type  of  registration  are  a  few  pixels. 


B  ©  0  \  If4  Screli  ra.058741 


Figure  17.  Image  to  image  registration  example  performed  with  a  few  ground  reference 

points  along  shoreline 
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After  each  of  the  first  three  images  were  co-registered  to  image  eight,  principal 
component  transforms  were  performed  on  each  image.  The  WorldView-2  multispectral 
images  consist  of  eight  distinct  spectral  bands;  Coastal  Blue,  Blue,  Green,  Yellow,  Red, 
Red-Edge,  NIRl,  and  NIR2  (DigitalGlobe,  2010b).  Eight  band  multispectral  images 
provide  increased  degrees  of  freedom  when  performing  principal  component  transforms 
of  the  data.  Figure  18  shows  each  of  the  eight  original  spectral  bands  for  the  first  image 
of  an  area  with  many  visible  waves.  Figure  19  shows  each  of  the  eight  principal 
component  bands  for  the  first  image.  Principal  component  four  highlights  the  waves  with 
distinct  narrow  lines  at  the  wave  crests  without  the  additional  detail  and  washed  out 
features  shown  in  the  other  principal  components. 


Figure  18.  Original  spectral  bands  one  through  eight  example 
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Figure  19.  Principal  component  bands  one  through  eight  example 

Taking  a  transect  of  the  inverse  of  the  principal  component  four  image  shows 
high  values  corresponding  to  wave  crest  position.  Figure  20  shows  a  transect  of  the 
principal  component  four  image  with  transect  profile. 


Figure  20.  Principal  component  four  inverse  image  with  transect  perpendicular  to  shore 

and  profde 
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Once  the  principal  components  transform  has  been  performed  on  each  image,  the 
waves  in  each  image  can  be  compared  to  each  other.  A  new  image  can  be  displayed  to 
show  a  change  detection  from  the  first  image  to  the  second  image.  Principal  component 
four  of  the  second  image  is  placed  into  both  the  red  and  green  inputs  for  an  RGB  triplet 
representation  of  the  data.  Principal  component  four  of  the  first  image  is  placed  in  the 
blue  input. 

Displaying  the  two  images  in  this  way  allows  visualization  of  the  position  of  the 
visible  waves  in  the  two  images  as  they  change  from  the  first  image  to  the  second  image. 
Figure  21  shows  the  resultant  change  detection  image.  The  blue  lines  represent  the 
highlighted  wave  crests  from  principal  component  four  in  the  first  image.  The  yellow 
lines  (the  additive  color  combination  of  red  and  green)  represent  the  highlighted  wave 
crests  from  principal  component  four  in  the  second  image.  The  change  detection  image 
represents  how  the  waves  have  moved  from  the  first  image  to  the  second  image  in  the 
time  between  the  two  images.  The  waves  have  moved  from  the  blue  line  positions  to  the 
yellow  line  positions  in  time  10.80  seconds.  This  allows  us  to  determine  the  change  in 
wave  position  from  one  image  to  the  next  image  in  the  time  between  images  which  gives 
us  the  wave  celerity.  The  distance  the  wave  moves  from  one  image  to  the  next  was 
determined  using  the  measurement  tool  in  ENVI. 
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Figure  2 1 .  Change  from  first  to  second  image  where  wave  moves  from  blue  to  yellow 

toward  the  shore 

We  also  looked  at  ways  of  highlighting  several  different  scene  components  using 
the  eight  multispectral  bands.  Figure  22  shows  part  of  the  scene  from  image  P002  row  1 
column  1.  The  image  was  displayed  with  bands  six,  five,  and  three  in  an  RGB  triplet 
representation  of  the  data.  Several  of  the  main  components  of  the  scene  are  highlighted 
including  from  left  to  right;  waves,  kelp  (red),  an  airplane  condensation  trail,  outflowing 
river  sediment,  and  the  shoreline  (white). 
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Worldview  2 
Camp  Pendleton 

2010  March  24-  190332Z 
P002  Row  1  Column  1 
RGB:  Bands  653 

k;  3367:5454 
36735151 


Figure  22.  Illustration  of  Camp  Pendleton  scene  components  using  multispectral  data. 

Components  of  the  scene  highlighted  include  from  left  to  right;  waves,  kelp  (red), 
an  airplane  condensation  trail,  outflowing  river  sediment,  and  the  shoreline 

(white) 


Figure  23  shows  the  same  scene  using  both  multispectral  and  panchromatic  image 
data.  Principal  component  zero  (intensity)  from  the  multispectral  image  of  P002  row  1 
column  1  was  compared  to  the  panchromatic  image  of  the  same  scene.  In  this  image,  the 
kelp  (red)  is  now  much  better  resolved  because  it  is  being  represented  with  the 
panchromatic  data  which  has  approximately  four  times  higher  resolution.  The  principal 
component  zero  data  shows  the  same  features  as  in  the  previous  figure,  but  the  features 
are  represented  with  different  colors.  Figures  22  and  23  are  illustrations  of  our  ability  to 
separate  scene  components  with  the  available  image  data. 
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Figure  23.  Illustration  of  Camp  Pendleton  scene  components  using  multispectral  and 

panchromatic  data 
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B.  WORLDVIEW-2  COLLECTION  SIMULATION 

A  simulation  of  the  Camp  Pendleton  area  image  collection  was  created  using 
Satellite  Tool  Kit  (from  AGI).  Figure  24  shows  a  movie  of  the  collection  simulation: 


Figure  24.  Satellite  Toolkit  simulation  of  WorldView-2  Camp  Pendleton  area  collection 
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V.  ANALYSIS 


A.  BATHYMETRY  FROM  IMAGERY 

In  order  to  determine  depth  relatively  close  to  the  shoreline  using  linear  dispersion 
theory,  we  used  the  first  three  images  in  the  data  set.  Each  of  the  first  three  images 
(P002,  P007,  and  P009)  were  registered  to  one  of  the  middle  images  (POOS).  Principal 
component  transforms  were  then  computed  for  each  co-registered  image.  Statistics  were 
also  calculated.  Table  5  shows  the  computed  statistics  for  the  first  image  P002: 

Table  5.  Statistics  of  the  first  image  computed  during  principal  components 

transformation 


Filename:  /Users/blmccart/Documents/ENVI  Data/10MAR24190332-M2AS_R1C1-052365928020_01_P002.TIF 

Dims:  Full  Scene  (75,656,250  points) 

Basic  Stats 

Min 

Max 

Mean 

Stdev 

Num 

Eigenvalue 

Band  1 

0 

1300 

356.05 

80.90 

1 

168,625.55 

Band  2 

0 

1499 

300.95 

75.65 

2 

27,243.87 

Band  3 

0 

2047 

362.68 

104.54 

3 

2,341.52 

Band  4 

0 

2035 

277.79 

97.06 

4 

222.48 

Band  5 

0 

1955 

185.21 

73.39 

5 

101.52 

Band  6 

0 

2047 

396.04 

185.39 

6 

89.68 

Band  7 

0 

2047 

454.60 

256.65 

7 

29.89 

Band  8 

0 

2047 

431.48 

245.80 

8 

7.85 

Covariance 

Band  1 

Band  2 

Band  3 

Band  4 

Band  5 

Band  6 

Band  7 

Band  8 

Band  1 

6,545.4 

5,988.3 

7,455.9 

6,058.4 

4,286.0 

5,577.2 

5,007.8 

4,641.0 

Band  2 

5,988.3 

5,722.3 

7,437.6 

6,260.3 

4,548.6 

5,412.7 

4,583.0 

4,205.3 

Band  3 

7,455.9 

7,437.6 

10,928.0 

9,814.8 

7,213.5 

11,953.8 

12,587.8 

11,811.0 

Band  4 

6,058.4 

6,260.3 

9,814.8 

9,421.3 

7,049.1 

12,194.9 

13,013.0 

12,328.0 

Band  5 

4,286.0 

4,548.6 

7,213.5 

7,049.1 

5,386.2 

8,646.5 

8,936.4 

8,463.2 

Band  6 

5,577.2 

5,412.7 

11,953.8 

12,194.9 

8,646.5 

34,370.4 

46,489.9 

44,520.5 

Band  7 

5,007.8 

4,583.0 

12,587.8 

13,013.0 

8,936.4 

46,489.9 

65,871.6 

62,956.2 

Band  8 

4,641.0 

4,205.3 

11,811.0 

12,328.0 

8,463.2 

44,520.5 

62,956.2 

60,417.1 

Correlation 

Band  1 

Band  2 

Band  3 

Band  4 

Band  5 

Band  6 

Band  7 

Band  8 

Band  1 

1 

0.98 

0.88 

0.77 

0.72 

0.37 

0.24 

0.23 

Band  2 

0.98 

1 

0.94 

0.85 

0.82 

0.39 

0.24 

0.23 

Band  3 

0.88 

0.94 

1 

0.97 

0.94 

0.62 

0.47 

0.46 

Band  4 

0.77 

0.85 

0.97 

1 

0.99 

0.68 

0.52 

0.52 

Band  5 

0.72 

0.82 

0.94 

0.99 

1 

0.64 

0.47 

0.47 

Band  6 

0.37 

0.39 

0.62 

0.68 

0.64 

1 

0.98 

0.98 

Band  7 

0.24 

0.24 

0.47 

0.52 

0.47 

0.98 

1 

1.00 

Band  8 

0.23 

0.23 

0.46 

0.52 

0.47 

0.98 

1.00 

1 

Eigenvector 

Band  1 

Band  2 

Band  3 

Band  4 

Band  5 

Band  6 

Band  7 

Band  8 

PC  Band  1 

-0.069 

-0.066 

-0.148 

-0.150 

-0.105 

-0.450 

-0.619 

-0.592 

PC  Band  2 

0.417 

0.416 

0.511 

0.436 

0.330 

0.053 

-0.209 

-0.214 

PC  Band  3 

0.645 

0.336 

-0.014 

-0.431 

-0.467 

-0.202 

0.131 

0.099 

PC  Band  4 

-0.362 

0.075 

0.522 

-0.085 

-0.067 

-0.506 

0.522 

-0.225 

PC  Band  5 

0.206 

0.042 

-0.465 

0.179 

0.520 

-0.599 

0.202 

0.194 

PC  Band  6 

-0.185 

0.020 

0.335 

-0.030 

-0.067 

-0.324 

-0.487 

0.711 

PC  Band  7 

-0.130 

0.275 

0.068 

-0.724 

0.583 

0.181 

-0.071 

-0.016 

PC  Band  8 

0.426 

-0.791 

0.332 

-0.192 

0.207 

-0.045 

0.019 

-0.020 

The  statistics  show  that  bands  7  and  8  have  very  high  covariance  compared  to  the 
other  covariance  values.  Figure  25  shows  a  plot  of  spectral  band  covariance  values. 
Band  to  band  correlation  was  also  computed.  Figure  26  shows  a  plot  of  spectral  band 
correlation  values. 
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Figure  25.  Covariance  for  the  spectral  bands  of  the  first  image 
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Figure  26.  Correlation  for  the  speetral  bands  of  the  first  image 
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Eigenvector  representations  of  each  of  the  eight  principal  component  bands  are 
shown  in  Figures  27  and  28.  Because  the  fourth  principal  component  is  best  for 
highlighting  wave  crests  and  giving  us  wave  position  at  a  specific  time,  we  were  most 
interested  in  its  eigenvector  representation.  Principal  component  four  is  highlighted  in 
Figure  27.  Principal  component  band  4  is  composed  of  large  parts  of  spectral  bands  three 
six  and  seven  (Green,  Red-Edge,  and  NIRl). 
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Figure  27.  Eigenvector  representation  of  principal  component  bands  one  through  four 

with  PC  Band  4  highlighted 
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Figure  28.  Eigenvector  representation  of  principal  component  bands  five  through  eight 

Once  principal  component  transforms  had  been  computed  for  each  co-registered 
image,  they  could  be  compared.  A  change  detection  image  was  created  for  the  first  two 
images.  The  change  detection  image  was  an  RGB  color  composite  with  principal 
component  four  for  the  second  image  in  the  red  and  green  components,  and  principal 
component  four  for  the  first  image  in  the  blue  component. 

In  the  change  detection  image,  blue  lines  show  the  position  of  the  wave  in  the  first 
image,  while  yellow  lines  show  the  position  of  the  waves  in  the  second  image.  The 
waves  move  from  their  positions  when  they  are  blue  to  their  positions  when  they  are 
yellow  in  the  time  between  the  two  images  (10.80  seconds).  Figure  29  shows  a  chip  of 
this  image.  In  the  second  change  detection  image,  the  waves  move  from  their  positions 
when  they  are  blue  to  their  positions  when  they  are  yellow  in  the  time  between  the  two 

images  (also  10.80  seconds).  Figure  30  shows  a  chip  of  this  image. 
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Figure  29.  First  image  to  second  image  change  detection 


Figure  30.  Second  image  to  third  image  change  detection 
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We  were  able  to  compare  the  two  change  detection  images  to  each  other  as  well. 
Another  change  detection  image  was  created  as  an  RGB  color  composite.  The  second  of 
the  original  change  detection  images  was  used  for  red  and  green.  The  first  of  the  original 
change  detection  images  was  used  for  blue.  In  this  new  change  detection  image  waves 
move  from  cyan  to  magenta  to  yellow  which  are  additive  color  combinations  of  the 
original  red,  green,  and  blue  inputs. 

Because  the  time  delta  between  first  and  second  images  is  the  same  as  the  time 
delta  between  second  and  third  images,  the  lines  lie  directly  on  top  of  each  other  and  the 
red,  green,  and  blue  inputs  are  all  color  mixed.  That  is  why  there  are  no  red,  green,  or 
blue  colors  in  the  resultant  image.  Figure  31  shows  a  chip  from  the  image.  From  the 
figure  you  can  see  that  there  is  an  additional  yellow  line  between  cyan  and  magenta. 
There  is  also  an  additional  cyan  line  between  magenta  and  yellow.  These  additional  lines 
are  caused  by  multiple  waves  overlapping.  Multiple  waves  overlapping  will  cause  more 
and  more  additional  lines  to  appear  as  more  data  sets  are  combined  in  this  way. 
Combining  more  than  two  sets  of  images  in  this  way  was  not  particularly  useful. 


Figure  3 1 .  First  to  second  to  third  image  change  detection  showing  only  additive  color 

mixed  wave  lines 
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The  measurement  tool  in  ENVI  was  used  to  measure  the  distance  the  wave 
travels.  The  depth  of  the  water  at  the  midpoint  between  the  blue  and  yellow  lines  in  the 
change  detection  image  was  determined  using  Equation  1  (Equation  4  can  also  be  used). 
Water  depth  was  determined  in  this  way  at  points  all  along  the  shoreline  in  the  change 
detection  images. 

Next,  it  was  necessary  to  compare  these  computed  depths  to  reference 
bathymetry.  Maps  were  created  of  the  Camp  Pendleton  area  using  ArcGIS  software.  The 
maps  allowed  accurate  alignment  of  image  data  with  reference  bathymetry.  The 
bathymetry  used  was  created  by  ESRI  and  modified  by  the  California  Department  of  Fish 
and  Game.  The  bathymetry  consists  of  ten  meter  contour  lines  to  six  hundred  meters 
depth.  Figure  32  shows  a  chip  from  a  change  detection  image  placed  in  a  map  with 
bathymetry  contour  lines  and  sediment  types. 


Figure  32.  Overlay  of  surveyed  bathymetry  to  depth  measured  change  detection  image 

This  process  works  well  for  the  waves  close  to  the  shore  (<1000  m  from  the 
shore)  where  principal  component  four  is  able  to  highlight  the  wave  crests.  For  waves 
farther  from  shore  (1000  m  to  3000  m),  principal  component  four  is  not  able  to  highlight 
the  wave  crests.  Another  process  must  be  used  in  order  to  determine  depth  farther  from 
shore.  The  wave  period  method  as  described  by  Williams,  1947  was  used  for  these 
waves. 
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Depth  curves  were  created  as  a  function  of  wave  period  and  wavelength  in  order 
to  perform  the  wave  period  determination  of  water  depth.  These  were  very  similar  to 
those  produced  by  Williams,  1947  (see  Figure  2).  Because  we  were  working  in  units  of 
meters,  a  set  of  depth  curves  was  created  using  units  of  meters  and  seconds.  The  figure 
was  created  to  allow  determination  of  depth  from  a  single  multispectral  image  of  the 
Camp  Pendleton  area.  Figure  33  shows  a  plot  of  depth  curves  as  a  function  of  period  and 
wavelength.  Wavelengths  for  waves  occurring  over  known  depth  postings  were 
measured  using  the  ENVI  measurement  tool.  Surface  gravity  waves  maintain  the  same 
period  as  they  propagate  but  their  celerity  and  wavelength  can  change.  Other 
wavelengths  can  then  be  measured  and  the  same  period  assumed.  The  depth  below  other 
waves  can  be  determined  by  which  depth  curve  intersects  the  period  and  wavelength  of 
the  wave  using  Figure  33. 
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Figure  33.  Curves  relating  wave  period  to  depth  and  wavelength  (wavelength  units  in 

meters) 


The  multispectral  image  data  were  represented  in  several  different  ways  to  see 
which  representation  could  best  highlight  waves  in  a  given  image  chip.  Figure  34  shows 
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one  representation  using  bands  seven,  six,  and  five  in  an  RGB  color  representation  of  the 
data.  These  bands  correspond  to  NIRl,  Red-Edge,  and  Red  bands.  These  bands  do  not 
penetrate  the  water,  and  this  results  in  increased  surface  wave  contrast  (Abileah,  2006). 
Different  band  combinations  and  image  enhancements  were  used  depending  on  the 
particular  image  chip.  Typically,  image  equalization  enhancement  did  the  best  job  of 
highlighting  the  waves  in  a  given  image  chip. 


Figure  34.  Deeper  water  waves  resolved  using  765  bands  in  RGB  triplet  representation 

We  now  had  two  different  methods  of  determining  water  depth  from  the  Camp 
Pendleton  area  imagery.  The  first  was  to  use  co-registered  sequential  images  transformed 
into  their  principal  components.  A  single  change  detection  image  could  then  be  created 
that  showed  the  change  in  principal  component  four  from  the  first  image  to  the  second. 
The  depth  below  the  halfway  point  between  the  two  wave  positions  could  be  computed 
using  the  wave  velocity  method  with  the  linear  dispersion  relation  equation.  This  method 
works  well  for  waves  having  a  distance  from  shore  of  less  than  a  thousand  meters  or  so. 

The  second  method  uses  a  single  multispectral  image  represented  as  an  RGB 
color  composite  of  three  of  the  higher  bands  and  equalization  enhancement.  The 
wavelength  above  a  known  depth  can  be  used  to  determine  wave  period.  Following  the 
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same  wave  direction  to  shore  we  can  assume  the  same  period  for  other  waves  of  differing 
wavelength.  Using  Figure  33,  we  can  then  determine  depth  for  the  other  waves  by 
measuring  their  wavelengths  and  assuming  the  same  wave  period. 

B.  COMPARISON  TO  REFERENCE  BATHYMETRY 

The  next  step  was  to  compare  these  calculated  depths  to  our  reference 
bathymetry.  The  WorldView-2  images  used  for  this  study  were  taken  between  12:03  PM 
and  12:06  PM  PDT.  The  National  Oceanic  and  Atmospheric  Administration  operates 
coastal  stations  where  tide  levels  are  verified.  The  closest  station  to  the  Camp  Pendleton 
area  is  in  La  Jolla  around  thirty  miles  south.  The  station  identification  number  is 
9410230.  The  NOAA  website  allowed  determination  of  the  verified  water  level  on  the 
day  and  time  of  image  acquisition.  Figure  35  shows  a  plot  of  water  level  versus  time  for 
March  24,  2010  at  La  Jolla,  California. 


Figure  35.  Tide  height  for  March  24,  2010  at  closest  station  to  Camp  Pendleton  (From 
[National  Oceanic  and  Atmospheric  Administration,  2010,  August  25]) 

Figure  35  shows  that  when  the  Camp  Pendleton  area  images  were  acquired  (12:03 
to  12:06  PM  PDT),  the  water  level  was  very  close  to  the  mean  low  water  level.  Our 
reference  bathymetry  is  also  relative  to  mean  low  water  level  so  there  is  no  real  water 
level  offset  that  needed  to  be  done  to  compare  calculated  depth  to  the  reference 
bathymetry. 
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The  first  change  detection  image  was  carefully  combined  with  the  map  layer 
containing  the  reference  bathymetry.  A  perpendicular  line  was  drawn  on  the  image  from 
the  shoreline  out  into  the  water  where  it  would  cross  as  many  contour  lines  as  possible. 
Figure  36  shows  the  resultant  image  with  perpendicular  drawn  to  the  shoreline.  The 
ENVI  measurement  tool  was  used  to  measure  from  the  shoreline  to  where  each  contour 
line  crossed  the  perpendicular  line.  This  produced  a  reference  bathymetry  depth  profile 
curve. 


Figure  36.  Bathymetric  contours  with  perpendicular  drawn  to  shoreline  overlayed  on  top 

of  change  detection  image 

Once  the  reference  bathymetry  depth  profile  was  determined,  we  determined 
depth  for  other  waves  along  the  perpendicular  line.  For  waves  less  than  about  five 
hundred  meters  from  the  shoreline,  the  change  detection  image  was  used  to  determine 
depth  using  the  wave  velocity  method.  Only  four  points  along  the  perpendicular  line 
were  measured  since  only  four  were  easily  distinguished  in  the  image. 

A  multispectral  image  representation  of  the  data  was  created  in  a  new  window 

and  the  two  images  were  linked  together  using  ENVI.  In  this  way  we  were  able  to  use 

the  reference  image  containing  the  reference  bathymetry  and  perpendicular  line  with  a 

single  image  of  the  scene.  The  wavelength  of  a  wave  crossing  the  twenty  meter  depth 
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contour  was  measured  and  its  period  was  determined  from  Figure  33.  The  wavelength  of 
waves  along  the  perpendicular  line  toward  the  shore  were  then  measured  and  using  the 
same  period  as  the  wave  at  twenty  meters  depth,  the  depths  of  the  midpoint  between  the 
waves  crests  was  determined.  For  each  depth  point,  the  distance  to  shore  was  determined 
using  the  ENVI  measurement  tool.  Figure  37  shows  the  resultant  depth  profile. 


Camp  Pendleton  Nearshore  Depth  Profile 
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Figure  37.  Nearshore  depth  profile  with  data  points  from  wave  velocity  and  wave  period 

methods 

The  red  triangles  represent  the  data  points  determined  through  the  use  of  the 
change  detection  image  and  the  wave  velocity  method  using  the  linear  finite  depth 
dispersion  relation  for  surface  gravity  waves.  The  green  circles  represent  the  data  points 
determined  through  the  use  of  the  wave  period  method  using  the  linear  dispersion 
relation.  From  about  six  hundred  meters  to  twelve  hundred  meters  from  the  shoreline, 
there  is  an  airplane  condensation  trail  which  obscured  the  waves  and  we  were  not  able  to 
get  depth  data  points  in  that  area.  For  distances  greater  than  about  thirty  five  hundred 
meters  from  shore  waves  could  not  be  very  well  resolved  so  we  were  not  able  to  get 
depth  data  points  beyond  that  point. 
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Using  the  linear  dispersion  theory  has  a  maximum  predicted  error  of  42% 
associated  with  not  taking  into  account  the  non  linearities  (primarily  the  height  of  the 
waves)  associated  with  surface  gravity  wave  dispersion  theory  (Stockdon  &  Holman, 
2000).  The  maximum  predicted  error  results  in  an  overprediction  of  water  depth.  Figure 
38  shows  the  depth  profde  plot  further  zoomed  in  on  the  data  points.  The  data  points 
associated  with  the  wave  velocity  method  of  depth  determination  have  +42%  error  bars 
on  them.  The  data  points  associated  with  the  wave  period  method  of  depth  determination 
have  ±  15%  error  bars  on  them  since  they  are  farther  from  shore  where  the  accuracy  of 
the  reference  bathymetry  should  be  higher.  These  error  bars  were  used  since  they 
allowed  most  of  the  points  greater  than  six  hundred  meters  from  the  shore  to  intersect 
with  the  reference  bathymetry  line. 

Depth  accuracy  of  the  data  points  is  not  readily  apparent  since  the  reference 
bathymetry  has  been  linearly  interpolated  between  the  few  data  points  available.  Given 
the  sparse  reference  bathymetry  it  is  possible  that  our  calculated  depths  are  more  accurate 
than  the  reference  bathymetry.  The  bottom  sediment  composition  in  the  nearshore  area 
of  this  location  consists  primarily  of  fine  sand  and  silt.  Seasonal  changes  in  ocean 
currents  and  tide  levels,  as  well  as  severe  storms  could  cause  changes  to  depth  because  of 
the  variable  nature  of  the  bottom  sediments.  It  is  possible  that  this  could  cause 
discrepancies  between  the  reference  bathymetry  and  the  measured  depths,  especially  in 
the  surf  zone.  High  resolution  bathymetry  (like  high  resolution  imagery)  necessarily 
equates  to  lower  area  coverage.  In  order  to  get  higher  resolution  bathymetry  covering 
the  entire  area  in  our  scene,  we  would  need  to  perform  our  own  sonar  bathymetric  survey. 
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Camp  Pendleton  Nearshore  Depth  Profile 


Distance  From  Shore  [m] 

Figure  38.  Nearshore  depth  profile  with  error  bars  (plus  42%  for  wave  velocity,  and 

plus/minus  15%  for  wave  period) 
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VI.  SUMMARY 


WorldView-2  multispectral  imagery  of  the  coastal  area  near  Camp  Pendleton, 
California  were  used  to  determine  ocean  depth  using  the  linear  finite  depth  dispersion 
relation  for  surface  gravity  waves.  The  high  spatial  resolution  (fifty  centimeters) 
available  from  the  panchromatic  imagery  was  not  necessary  for  resolving  waves  since 
waves  are  much  larger  than  this.  The  increased  spectral  resolution  from  the  eight  bands 
composing  the  multispectral  imagery  provided  more  degrees  of  freedom  which  could  be 
used  to  perform  the  measurements. 

The  multispectral  images  were  co-registered  and  principal  component  transforms 
were  performed  on  them.  Principal  component  four  was  found  to  very  clearly  highlight 
wave  crests  for  waves  within  one  thousand  meters  from  shore.  Creating  a  change 
detection  image  with  principal  component  four  for  the  second  image  represented  by  red 
and  green,  and  principal  component  four  for  the  first  image  represented  by  blue  in  an 
RGB  triplet  representation  of  the  data  produced  a  single  image  containing  both  spatial 
and  temporal  information.  In  the  change  detection  image,  waves  moved  from  blue  to 
yellow  in  the  time  between  images.  This  allowed  us  to  accurately  measure  the  distance 
traveled  by  the  wave  in  the  time  between  images,  which  gave  us  wave  celerity.  Depth 
inversion  of  nearshore  waves  could  be  performed  once  wave  celerity  was  known. 

For  waves  between  one  thousand  and  thirty-five  hundred  meters  from  shore, 
another  method  for  determining  depth  was  needed  since  principal  component  four  no 
longer  highlighted  wave  crests  in  this  region.  The  multispectral  imagery  was  represented 
as  an  RGB  composite  image  with  longer  wavelength  bands  and  was  equalization 
enhanced.  This  imagery  allowed  us  to  use  the  wave  period  method  to  determine  ocean 
depth.  The  calculated  depth  data  points  were  compared  to  reference  bathymetry.  Error 
bars  of  -1-42%  were  used  on  points  determined  with  the  wave  velocity  method.  Error  bars 
of  ±15%  were  used  on  points  determined  with  the  wave  period  method.  Given  the  sparse 
reference  bathymetry  it  is  possible  that  our  calculated  depths  are  more  accurate  than  the 
reference  bathymetry,  especially  in  the  surf  zone. 
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VII.  CONCLUSIONS  AND  FUTURE  WORK 


Nearshore  ocean  depth  was  successfully  determined  using  multispectral  imagery 
acquired  by  the  WorldView-2  satellite  of  the  coastal  area  near  Camp  Pendleton, 
California.  Image  processing  techniques  were  used  to  enhance  wave  position  and  allow 
measurement  of  wavelengths.  Depth  inversion  was  performed  using  both  the  wave 
velocity  method  and  the  wave  period  method.  The  computed  depth  data  points  compare 
favorably  to  reference  bathymetry. 

This  work  shows  that  it  is  feasible  to  use  space-based  multispectral  sensors  to 
determine  nearshore  ocean  depth  in  denied  waters  where  recent  sonar  bathymetric  survey 
data  are  not  available.  The  linear  finite  depth  dispersion  relation  for  surface  gravity 
waves  requires  only  that  the  top  of  the  water  waves  can  be  resolved  as  opposed  to  the 
transparency  method  of  depth  determination. 

Future  work  on  this  subject  should  involve  an  extensive  high  resolution  sonar 
bathymetric  survey  of  the  coastal  region  in  the  scene  of  interest.  This  will  allow  better 
determination  of  computed  depth  accuracy  when  compared  to  recent  high  resolution 
bathymetry  data. 

The  linear  dispersion  relation  is  capable  of  providing  fairly  accurate  depth 
determination  from  space  based  imagery.  The  use  of  the  non-linearity  associated  with 
wave  height  would  increase  the  accuracy  of  depth  inversion,  but  the  wave  trough  to  crest 
distance  is  very  small  compared  to  the  wave  to  spacecraft  distance.  This  may  make  it 
difficult  to  determine  wave  height  from  space. 

The  time  interval  between  images  should  be  closer  to  five  seconds  in  order  to 
increase  accuracy  without  adversely  affecting  wave  resolution.  Because  even  a  small 
amount  of  cloud  cover  will  affect  the  ability  to  resolve  waves,  more  images  (twenty  to 
thirty),  should  be  taken  of  the  same  scene  at  five  second  intervals.  Exhaustive  depth 
computations  using  many  images  of  the  scene  can  be  used  to  create  a  depth  data  point 
cloud  along  a  depth  profile  line  common  to  all  images.  This  should  help  to  establish  the 
accuracy  of  this  method. 
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